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Abstract. In cold dark matter cosmological models, the first stars to form are believed to do so within small protogalaxies. We 
wish to understand how the evolution of these early protogalaxies changes once the gas forming them has been enriched with 
small quantities of heavy elements, which are produced and dispersed into the intergalactic medium by the first supemovae. 
Our initial conditions represent protogalaxies forming within a fossil H II region, a previously ionized region that has not 
yet had time to cool and recombine. We study the influence of low levels of metal enrichment on the cooling and collapse 
of ionized gas in small protogalactic halos using three-dimensional, smoothed particle hydrodynamics (SPH) simulations 
that incorporate the effects of the appropriate chemical and thermal processes. Our previous simulations demonstrated that 
for metallicities Z < 10^"' Zq, metal line cooling alters the density and temperature evolution of the gas by less than 1% 
compared to the metal-free case at densities below Icm^"' and temperatures above 2000 K. Here, we present the results of 
high-resolution simulations using particle splitting to improve resolution in regions of interest. These simulations allow us to 
address the question of whether there is a critical metallicity above which fine structure cooling from metals allows efficient 
fragmentation to occur, producing an initial mass function (IMF) resembling the local Salpeter IMF, rather than only high- 
mass stars. 
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INTRODUCTION 

Population III stars are the first potential producers of 
UV photons that can contribute to the reionization pro- 
cess and are the first producers of the metals required for 
the formation of population II stars. Metals produced by 
the first stars will also be injected into some fraction of 
the ionized volume, and so will be present in gas falling 
into new or existing protogalactic halos. The question 
then arises as to how this low level of metal enrichment 
affects the ability of the gas to cool and collapse. 

In [1], the collapse of cold, metal-enriched gas in a 
top-hat potential was simulated. These simulations in- 
cluded the effects of atomic fine structure cooling, but did 
not include cooling from H2. In the absence of molecu- 
lar cooling, fragmentation suggestive of a modern IMF 
was found to occur only for metallicities above a thresh- 
old value of Z ~ lO^^'^Zp^. However, the authors noted 
that the neglect of molecular cooling could be signifi- 
cant. In |2], Omukai et al. argued, based on the results of 
their detailed one-zone models, that molecular cooling 
would indeed dominate the cooling over many orders of 
magnitude in density. Also, by starting with cold gas, the 
authors of 1 1] implicitly assume that no extra entropy or 
energy has been added to the gas during its enrichment. 



although this is unlikely to be the case |3]. 

It is therefore important to consider whether a metal- 
licity threshold appears in simulations with different ini- 
tial conditions. In this paper, we examine a different, still 
idealized, set of initial conditions, using a more detailed 
treatment of the cooling and chemistry of the gas, and 
find no threshold. We argue that the transition from the 
primordial to the modern IMF therefore depends on the 
conditions under which stars form as much as on the 
metal abundances present. 

NUMERICAL APPROACH 

To help us to assess the influence of metals on the frag- 
mentation properties of gas in small protogalactic halos, 
we have performed a number of numerical simulations. 
Here we present a portion of our results. Detailed dis- 
cussion of the full set of simulations can be found in 
1(4)] (hereafter paper V). Our study requires us to follow 
the collapse of the gas over many orders of magnitude 
in density, and therefore the problem is well-suited for 
the use of a Lagrangian numerical method. We use SPH, 
a Lagrangian method described in isl], |@] and 101, in 
which the fluid is followed with an ensemble of particles. 



Fluid properties at each point are computed by averaging 
over neighboring particles. The calculations presented in 
this paper use version 1 of the parallel SPH code GAD- 
GET 1 8]. In order to achieve a higher mass resolution, 
we refine the mass of the gas particles using a method 
known as particle splitting |9]. In 1 10], hereafter paper II, 
we presented low-resolution simulations that show that 
the density increases most rapidly close to the center of 
the dark matter halo. On-the-fly splitting with two lev- 
els of refinement at two different radii provides us with 
the highest mass resolution at the region of interest. In 
the region of collapse we can resolve Jeans masses down 
to 1.5 Mq, as compared to Mres — 200 Mq in the low 
resolution simulations. We discuss the details of the im- 
plementation of the particle splitting algorithm in paper 
V. 

As we want to represent gas that has collapsed beyond 
the resolution limit of the simulation in a numerically 
robust manner, we have modified the code to allow it 
to create sink particles - massive, non-gaseous particles, 
designed to represent dense cores, that can accrete gas 
from their surroundings but otherwise interact only via 
gravity L 1 1,1 . Sink particles are created once the density 
rises above 10^ cm^^* and are endowed with an accretion 
radius of 0.3pc. On every time step, any gas within this 
accretion radius that is gravitationally bound to the sink 
particle is accreted by it. The design and implementation 
of our sink particle algorithm is discussed in more detail 
inLL2J. 

CHEMISTRY AND COOLING 

We have further modified GADGET to allow us to fol- 
low the non-equilibrium chemistry of the major coolants 
in both primordial and low-metallicity gas. Our chemi- 
cal model is presented in li 1 3i1 and LI 4,1 (hereafter papers 
I and III). Provided that carbon and oxygen are amongst 
the most abundant metals, the major coolants will be 
largely the same as in local atomic and molecular gas, 
namely H2, HD, C, C+, O, Si, Si+, CO, OH and H20II. 
We therefore follow the abundances of these ten species, 
together with an additional 29 species that play impor- 
tant roles in determining the abundances of one or more 
of these coolants. Our chemical network contains a to- 
tal of 189 collisional gas-phase reactions between these 
39 species, as well as 7 grain surface reactions, 43 reac- 
tions involving the photoionization or photodissociation 
of chemical species by ultraviolet radiation, and 8 reac- 
tions involving cosmic rays. For simplicity, in the sim- 
ulations presented in this paper we do not include the 
effects of dust, UV radiation or cosmic rays, and so use 
a simplified version of the model that contains only the 
collisional reactions. The thermal evolution of the gas in 
our simulations is modelled using a cooUng function that 



includes the effects of atomic fine structure cooling from 
C, C+, O, Si and Si+, rotational and vibrational cooling 
from H2, HD, CO and H2O, Lyman-a cooling, Comp- 
ton cooling, and H+ recombination cooling, as well as 
a number of other processes of lesser importance. For 
additional details of our chemical networks and cooling 
function, please consult papers I and III; further details 
of their implementation within Gadget are also given in 
paper II. 

INITIAL CONDITIONS 

Our initial conditions are based on those used in paper 
II, although collapse is followed to much higher density 
in the current work. We study protogalaxies forming 
from fully ionized gas with initial temperature T^, = 
10^ K, representing a fossil H II region |3] polluted 
by supernovae from the ionizing object. We model one 
such halo by using a fixed background potential with a 
spherically symmetric density profile |l5|]. We study the 
evolution both of zero-metallicity, primordial gas and of 
metal-enriched gas with a metallicity Z — lO^^Z©. To 
simplify the discussion of our simulation results, we take 
as a fiducial example a halo with a dark matter mass 
Mdm = 7.8 X 10^ Mq, an initial redshift z, = 25, and a 
spin parameter A = 0.05. For this example halo, the virial 
temperature T^^y — 1900 K, the virial radius ryh- = 0.1 kpc, 
and the truncation radius rt = 0.49 kpc, where both radii 
are given in physical units. The scale radius of this 
example halo is 29 pc, and the full computational volume 
is a box of side length Ikpc. The gas mass is Mg — 

0.19Mdm. 

RESULTS AND DISCUSSION 

In the simulations presented here we focus on gas with a 
metallicity Z = 10^^ Z©. This metallicity is at the upper 
limit of the range of values proposed for the so-called 
critical metallicity Z^nt, the value of the metallicity at 
which efficient fragmentation and low-mass star forma- 
tion is hypothesized to first occur 0, US [l3l- It is 
also comparable to the globally averaged metallicity pro- 
duced by the sources responsible for reionization [see 
e.g.lil. 

In Figure [T] we compare the temperature-density evo- 
lution of the run with primordial gas and the run with 
enriched gas. The two runs differ by less than 10% up to 
densities of n = 10^ cm^^. This result demonstrates that 
the cooling of the gas is barely influenced by the pres- 
ence of metals: fine-structure cooling contributes only 
marginally to the total cooling rate. H2 is the dominant 
coolant, rather than metal fine structure lines. At first 
sight, the fact that fine structure cooling from metals has 
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FIGURE 1. Gas temperature vs. number density for the 
runs with metallicities 0.0 {solid line), 10^^ (dotted line), and 
10^' Zq (dashed line). The spin parameter of the runs is 0.05. 
The time of the plot is half a Hubble time. The thin lines show 
the 1-<T deviation. 

little impact on the thermal or dynamical evolution of 
the gas at metallicities below 0.1 is somewhat sur- 
prising, given that it was previously found that gas with 
a metallicity of only Z = 10^^ Zr.., could cool rapidly and 
fragment even in the complete absence of molecular hy- 
drogen However, this difference in conclusions ap- 
pears to be a consequence of the different initial condi- 
tions used in the two sets of simulations, and the differ- 
ence in the physics included. The simulations presented 
in iQl] neglected molecular cooling, and started with cold 
(^initial = 200 K), low ionization (Xf, — 10^^) gas in a flat, 
top-hat dark matter potential lacking a central condensa- 
tion, although with local density perturbations. In these 
simulations, the gas cools, collapses and fragments only 
if the metallicity is high enough to render the cooling 
time shorter the the local dynamical timescale. With our 
arguably more realistic initial conditions, and, crucially, 
with our inclusion of H2, we find no metallicity threshold 
- H2 cooling alone is sufficient to allow the gas to col- 
lapse, but the collapsing gas does not fragment (see Fig- 
ure|2]and 1 19]). In |4], we present a set of additional runs 
with metallicities up to Z = lO^^Z© that confirm this re- 
sult. We conclude that the transition from the primordial 
to the modern-day IMF therefore depends on the condi- 
tions under which stars form as much as on the metal 
abundances present. The actual conditions that are ap- 
propriate for second generation star formation still need 
to be determined in detail by observation and modeling 
of galaxy formation. 
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FIGURE 2. Cut at z = 0.5 kpc through the gas number den- 
sity for the run with zero metallicity (left panel) and the run 
with a metallicity of Z = IO^-'Zq (right panel). We show a box 
with a size of 20 pc. 
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